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Q ' Abstract 

^ ■ A recent description of the motion of atoms in a classical monatomic 

system in liquid and supercooled liquid states divides the motion into 
^ , two parts: oscillations within a given many-particle potential valley, 

OO I and transit motion which carries the system from one many-particle 

valley to another. Building on this picture, we construct a model 
^SJ . for the trajectory of an average atom in the system. The trajectory 

I consists of oscillations at the normal-mode distribution of frequen- 

cies, representing motion within a fluctuating single-particle well, in- 
terspersed with position- and velocity-conserving transits to similar 
I adjacent wells. For the supercooled liquid in nondiffusing states, the 

model gives velocity and displacement autocorrelation functions which 

■ exactly match those found in the many-particle harmonic approxima- 
C , tion, and which are known to agree almost precisely with molecular dy- 
namics (MD) simulations of liquid Na. At higher temperatures, by al- 
lowing transits to proceed at a temperature-dependent rate, the model 
gives velocity autocorrelation functions which are also in remarkably 

^ I good agreement with MD simulations of Na at up to three times its 

■ melting temperature. Two independent processes in the model relax 
velocity autocorrelations: (a) dephasing due to the presence of many 
frequency components, which operates at all temperatures but which 
produces zero diffusion, and (b) the transit process, which increases 
with increasing temperature and which produces diffusion. Compared 
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to several treatments of velocity autocorrelations based on instanta- 
neous normal modes, the present model offers an advantage: It pro- 
vides a single-atom trajectory in real space and time, including transits 
and valid for arbitrary times, from which all single-atom correlation 
functions can be calculated, and they are also valid at all times. 

1 Introduction 

The motion of atoms in a liquid can be divided into two constituent parts: (a) 
oscillation in a valley of the liquid's many-body potential and (b) transits be- 
tween many-body valleys. The latter process is responsible for self- diffusion. 
As we showed in a previous paper [|T|, which drew upon earlier results of 
Clements and Wallace 0, |^, the former motion can be modeled very pre- 
cisely under the assumption that the valleys are nearly harmonic, with the 
majority of valleys (the random valleys) sharing a common spectrum of fre- 
quencies. Specifically, a purely harmonic model provides an extremely ac- 
curate formula for Z{t), the normalized velocity autocorrelation function, in 
the nondiffusing regime. Here we will use this picture and our previous work 
to justify a "mean atom trajectory" model, a single-atom model that approx- 
imates the behavior of an average atom in the liquid and correctly reproduces 
its nondiffusing behavior. Then we will introduce a simple intuitive account 
of the transit process that allows us to extend the model to the self-diffusing 
regime. In Sec. ^ we develop this model and explain how it is used to cal- 
culate Z{t) for a diffusing liquid. Then we fit the model to MD simulations 
at various temperatures in Sec. |^, and we comment on the quality of the 
results. In Sec. ^ we compare the present mean atom trajectory model with 
previous work based on Instantaneous Normal Modes (INM), and with an 
earlier independent atom model 0, and we summarize our conclusions. 

2 The Mean Atom Trajectory Model 
2.1 General comments 

To form an appropriate basis for our model, we must begin with some initial 
reasonable approximations about the nature of the valleys and transits in 
the real liquid. As mentioned in the Introduction, available evidence sug- 
gests that the valleys are nearly harmonic, and we will continue to assume 
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that here. Further, we will assume that transits between valleys occur in- 
stantaneously and are local in character; that is, each transit involves only a 
few neighboring atoms. 

Now for the nondiffusing supercooled liquid, there are no transits, and 
as we've seen the system's motion is accurately expressed in terms of the 
harmonic normal modes. But it is also legitimate to consider this motion 
from the point of view of a single atom, and when we do so we see that each 
atom moves along a complicated trajectory within a single-particle potential 
well which fluctuates because of the motion of neighboring atoms, but whose 
center is fixed in space. An important timescale for this motion is the single- 
atom mean vibrational period r, which we take as 27r/c(Jrms) where cUrms is 
the rms frequency of the set of normal modes. Let us now follow the single- 
atom description as the temperature is increased. Once the glass transition is 
passed, the atom will begin to make transits from one single-particle well to 
another 0, and even before the melting temperature is reached, the transit 
rate will be on the order of one per mean vibrational period. Since each atom 
has approximately ten neighbors, roughly ten transits will occur in its im- 
mediate vicinity every period, changing the set of normal mode eigenvectors 
each time; thus we conclude that a decomposition of the motion into normal 
modes will not be useful when the liquid is diffusing, so we have no choice 
but to follow a single-atom description for the liquid state Q]. That being 
the case, we shall start from the beginning with a single-atom description, 
in order to construct a unified model for diffusing and nondiffusing motion 
alike. 

2.2 Nondiffusing regime 

Our starting point will be in the nondiffusing regime, where a normal mode 
analysis is still valid. As shown in the ith coordinate of the i^'th atom in 
an iV-body harmonic valley can be written 



where the WKi,x form a 3A^ x 3A^ orthogonal matrix, the uj\ are the frequencies 
of the normal modes, and the ax are the amplitudes of the modes. Three 
of the modes have zero frequency and correspond to center of mass motion; 
here we demand that the center of mass is stationary, so the system has only 
3A^ — 3 independent degrees of freedom, the zero frequency modes are absent. 
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and the sum over A runs from 1 to 3A^ — 3. To make this an equation for a 
"mean" atom, we first drop the index K: 

i't) = J2 ^^aOa sin(u;At + ax). (2) 

A 

Now we must reinterpret the Wix since they no longer form an orthogonal (or 
even square) matrix. Let 

Wx = WixX + W2xy + WsxZ (3) 

so 

u{t) = sin(t^Ai^ + "a)- (4) 

A 

We will ultimately consider situations in which the well center is allowed to 
move, so let r{t) be the atom's position, R be the location of the center of 
the well, and u{t) be the atom's displacement from the well center; then 

r{t) = R + Ywxaxsin{ujxt + ax) (5) 

A 

with velocity 

v{t) = J2 i^xo-x^x cos{uJxt + ax). (6) 

A 

We now have the basic formula, but we must decide how to assign values 
to the Wx, ax, and ax. Let's do this by calculating Z{t), the velocity auto- 
correlation function, in this model and comparing it to the harmonic result 
derived in [|l|. 

Z{t) = l{v{t).v{0)) 

= ^ H ■ wyaxax'UJxUJx' {cos{uJxt + ax) cos{ay)) (7) 

^ AA' 

Let's assign the ax randomly and average over each ax separately; then if 

{cos{uJxt + ax) cos{ax')) = {cos{uJxt + oa)) (cos(a;A')) = (8) 
but if A = A', then 

(cos(u;At + ax) cos(aA')) = (cos(u;At + ax) cos(aA)) = 7: cos{uJxt). (9) 
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Thus the A 7^ A' terms in Eq. (|^ are ehminated and we have 

Z{t) = ^J2 I^Apa^cu^ cosM). (10) 

Note that the formula for Z{t) in the harmonic theory also lacks off-diagonal 
terms |jl|], but for a different reason: The orthogonality of the matrix WKi,\ 
removes them in that case. Now let us make the assignment 

1 . 

Wx = I W\ (11) 

where Wx is a randomly chosen unit vector; then 

^(^) = ^]^E«A^ACOs(^At). (12) 

This is the same expression for Z{t) that one derives in the harmonic model 
(see Eq. (10) of fl]]) for any distribution of normal mode amplitudes ax, not 
just the thermal equilibrium distribution. Thus our model with these choices 
of ax and Wx correctly reproduces Z{t) for any equilibrium or nonequilibrium 
ensemble in the harmonic theory. To recover the equilibrium result, we make 
the final substitution 



with the result 

1 kT 

Z{t) = — — ycos(cuAt), (14) 

^ ^ 3iV-3 M Y 

which is Eq. (13) from 0]. Thus our model assumes that the motion of a 
mean nondiffusing atom in thermal equilibrium at temperature T is given by 



1 2kT 

""^^^ " ^ ^ TW^ild ^ ''''^''"^ + ""^^ ^^^^ 

where the phases ax and unit vectors Wx are randomly chosen. 

By construction this model gets the same result for Z{t), and thus (f^), 
as the harmonic model; does it correctly reproduce any other functions? We 
can check by calculating the analog of Z{t) for positions, {u{t) ■ w(0)). Using 
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Eq. (|T]) for UKiit) and calculating for the harmonic model as in [Q, and using 
Eq. for u{t) and calculating as above, one finds in both cases that 



Ht)-n(o)) 



2N -2 



cos{uJxt). 



(16) 



Again the substitution ax = / Muj\ recovers the correct thermal equi- 

librium result. So this model correctly reproduces {u{t) ■ it(0)) and {v?) as 
well. Notice that because we have a closed form for u{t), all of the correla- 
tion functions calculated so far also have a closed form. Once we introduce 
transits, this will no longer be the case. 

While this model compares well with the harmonic model, one might 
wonder how well it compares with molecular dynamics (MD) results. In Fig. 
|l] a graph of f ^ as a function of t for the model in equilibrium at 6.69 K is 
compared to a randomly chosen particle from an MD run of liquid Na also 
at 6.69 K, a temperature at which it is known the sample is nondiffusing 
(see for details). Note that in both graphs v"^ has approximately the same 
amplitude, and the gaps between peaks are roughly the same size, indicating 
oscillations at frequencies in the same ranges. Thus not only functions of the 
motion but the motion itself shows strong qualitative agreement with MD in 
the nondiffusing regime. 

The MD system with which we will be comparing our model has = 500 
particles, so it has 1497 normal mode frequencies. (Remember that the zero 
frequency modes are removed at the start.) When we introduce transits, we 
will have to evaluate Z(t) numerically, and for speed of computation we would 
like to use only a representative subset of the normal mode frequencies, say 75 
instead of all 1497. The original set of frequencies is determined as described 
in We decided which subset to use by calculating three moments of the 
full frequency distribution defined below: 



CU2 



1/2 



exp(ln(c^A)) 



-1/2 



(17) 



Note that by this definition c<Jrms = \fs/5u!2- We then calculated the same 
three moments for several sets of 75 frequencies evenly spaced throughout 
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Figure 1: for a randomly chosen atom in an MD simulation of Na at 
6.69 K (top) is compared with for the mean atom trajectory model in 
the nondiffusing regime (bottom) over 5000 timesteps of duration 6t — 1.4 x 
10~^^ s. Notice that both fluctuations have roughly the same magnitude and 
oscillate at approximately the same frequencies. 
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the full set, and we chose the set that best fit the moments to use in further 
calculations. This process is somewhat subjective, because skewing the sam- 
ple in favor of lower frequencies improves the accuracy of a;_2 but reduces 
the accuracy of a;2, and the opposite is true if one skews in favor of high 
frequencies. The table below shows the values of the moments for the full set 
and the reduced set of 75 that we ultimately chose. The frequencies are in 
units of 6t~^, where 6t = 1.4 x 10~^^ s is the timestep of our MD simulations 

i- 

Full Reduced 

0.02826 0.02824 

cuo 0.01807 0.01803 

uj_2 0.02101 0.02111 

To use the reduced set we must rewrite our formulas for r{t) and v{t) slightly. 
Let A be the total number of frequencies; then A = 3A^ — 3 for the full set 
so r(t) and v(t) become 

[3 [2kT ^ 

r{t) = il+ y-y^^^A^r^sin(^At + aA) 

/3 [2kT ^ 
X\/^5I^ACOs(a;At + aA)- (18) 

This form is also correct for the reduced set of frequencies ux, so this is the 
form we will use. As another check on the accuracy of our results with only 
75 frequencies, we recalculated the original J^cos^uxt) expression for Z{t) 
using both the full and reduced sets; their disagreement is at most 0.01 out 
to time 1000 6t. However, the discrepancy grows beyond that time, as the 
expression using the reduced set begins to experience revivals. Hence we will 
not consider Z{t) beyond that point. 

2.3 Diffusing regime 

To incorporate diffusion into the mean atom trajectory model, we rely on 
Wallace's notion of a single-particle transit 0, a nearly instantaneous tran- 
sition from one well to another. As discussed in we expect transits to 
be governed not by thermal activation (having enough energy to escape a 
fixed well) but by correlations (neighbors must be positioned properly for a 
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low-potential path to open between wells). We implement this property by 
having transits occur at a temperature dependent rate i^{T), so in a small 
time interval At the probability of a single transit is uAt. 

We model the transit process itself by assuming it occurs instantaneously 
in the forward direction; from this we can determine the parameters R, Wx, 
and a\ appearing in Eq. (|1^) after the transit in terms of the same quantities 
before the transit. Since the process is instantaneous, both r{t) and v{t) are 
the same before and afterwards. Let R^'^^"^'^ and R^^^'^^ be the well centers 
before and after the transit and let and ti^^*'^'' be the corresponding 

displacements from the well centers. Then r^'^^°'^'^ = ^after jj-Qp^igg 

^before ^ ^before ^ j^aftcr ^ ^after_ ^-^g^j 

To transit forward, we assume the center of the new well lies along the line 
between the old well center and the atom, but it lies on the opposite side 
of the atom from the old well center an equal distance away. This implies 

^after ^ _u^elore ^ 

^before _|_ ^before j^after ^before (20) 

with the result 

j^after ^ j^before ^ 2v}"'^°"' . (21) 

This determines the new well center in terms of the coordinates before the 
transit. As for the unit vectors Wx, since they are randomly generated and 
play no role in calculating Z{t), we have decided to leave them unchanged 
by transits. Finally, we have the phases a\. We must use these to implement 
the relations 

after before „, after before (00\ 

Lb — LL , U — U I 

which we have assumed above. Since u{t) is a sum of sines and v{t) a sum 
of cosines, the simplest way to change the sign of u{t) while preserving that 
of v{t) is to reverse the signs of the arguments {uixt + ax) in Eq. ([18|). Let 
the transit occur at time to! then 

u^xto + af = - {ujxto + aY^n (23) 

so 

„ after o , j. „ before fr\A\ 

= -2uJxto - «A ■ (24) 
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Thus, in this model a transit is implemented at time to by leaving the W\ 
alone and making the substitutions 



R R + 2u{to) 

ax -^iJxto - ax- (25) 

This conserves r, reverses the sign of u, and conserves v. 

Now our mean atom trajectory model consists of nondifTusive motion 



between transits as given by Eq. (|T8D, with a given probability in each small 
time interval that R and the phases ax will be replaced with new values as 
determined in Eq. (^Sf). The addition of transits means that we no longer 
have closed form expressions for r{t) and v{t) for all times, so we have no 
closed form expression for Z{t); but this model can be implemented easily on 
a computer in a manner analogous to an MD calculation, and in that way we 
can calculate autocorrelation functions. We turn to that calculation next. 



2.4 Evaluating Z(t) in the diffusing regime 

To calculate Z{t), we select a value for the rate u, generate a random set of 
wx and ax, and use Eq. ( |T8| ) to calculate r(t) and v(t) from t = tot = tmax 
in increments of 6t, where the criterion for choosing tmax is discussed below 
and 6t is the timestep used in our MD simulations (defined in Subsection 
p.2|) . At each timestep, we check to see if a transit occurs, and if so we 
implement Eq. (|25|) and continue with the new R and ax- We then calculate 
Z{t) using the formula 

Zit) = _ E ^(^ + 1') ■ v{t') (26) 

■Jl^Jmax I) -\- o ^,_Q 

and normalizing. This equation is a modified form of the expression used to 
calculate Z{t) in MD; notice that the average over t' has the same effect as 
averaging separately over each phase ax that appears in the velocity vectors. 
Just as in MD, we want to average over a large data set, so we require 
^max >> t; we have chosen tmax = 20 million timesteps and we calculate Z{t) 
only to t = 1000 timesteps. We estimate the total error from using only a 
subset of all 1497 frequencies and the finite size of the data set to be at most 
0.01; in particular, when z/ = the calculation converges to the closed form 
result X) cos{uJxt) to this accuracy. 
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3 Comparison with MD 



The MD setup with which we compared the predictions of this model is 
the one described in = 500 atoms of Na move under the influence 

of a highly realistic pair potential with the timestep 5t = 1.4 x 10~^^ given 
in Subsection |2.2| . We performed equilibrium runs of the system at 216.3 K, 
309.7 K, 425.0 K, 664.7 K, and 1022.0 K, all temperatures at which the system 
is diffusing. Since Tm = 371.0 K for Na at this density, our simulations range 
from the supercooled regime to nearly three times the melting temperature. 
We then ran the model for various values of z/, adjusting until the model 
matched the value of the first minimum of Z{t) at each temperature. The 
values of u that we fit for all temperatures are given below, and the resulting 
Z{t) for each u is compared to the corresponding MD result in Figs. |^ through 
^ Here u is expressed in units of where r = 27r/a;rms is the single-atom 
mean vibrational period defined in Subsection p.l| . 

216.3 0.35018 
309.7 0.60276 
425.0 0.83985 
664.7 1.24858 
1022.0 1.68774 

Notice that in all cases u is of the same order of magnitude as r"^, indicating 
roughly one transit per mean vibrational period, as mentioned in Subsection 
and as predicted in p. 

The most obvious trend exhibited by Z{t) from the five MD runs is that 
its first minimum is rising with increasing T; as we mention below, this is 
the primary reason for the increasing diffusion coefficient D. Note that the 
model is able to reproduce this most important feature quite satisfactorily. 
In fact, all five fits of the model to the MD results capture their essential 
features, but we do see systematic trends in the discrepancies. First, note 
that the location of the first minimum barely changes at all in the model as 
I' is raised, but in MD the first minimum moves steadily to earlier times as 
the temperature rises. The first minimum occurs at a time roughly equal to 
half of the mean vibrational period (r = 287 6t in this system), so the steady 
drift backward suggests that the MD system is sampling a higher range of 
frequencies at higher T. Also, for the three lowest temperatures the model 
tends to overshoot the MD result in the vicinity of the first two maxima after 
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Figure 2: The model prediction for Z{t) at z/ = 0.35018 r ^ compared with 
the MD result for supercooled liquid Na at T = 216.3 K. 
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Figure 3: The model prediction for Z{t) at z/ = 0.60276 compared with 
the MD result for supercooled liquid Na at T = 309.7 K. 
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Figure 4: The model prediction for Z(t) at u 
the MD result for hquid Na at T = 425.0 K. 
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Figure 5: The model prediction for Z{t) at u 
the MD result for hquid Na at T = 664.7 K. 
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Figure 6: The model prediction for Z{t) aX u — 1.68774 r ^ compared with 
the MD result for hquid Na at T = 1022.0 K. 
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the origin, and at the highest two temperatures this overshoot is accompa- 
nied by a positive tail that is shghtly higher than the (still somewhat long) 
tail predicted by MD. These overshoots should clearly affect the diffusion 
coefficient D, which is the integral of Z{t). To check this, we calculated the 
reduced diffusion coefficient D, the integral of Z{t), which is related to D by 

D = %b. (27) 

The results are compared to the values of D calculated from the MD runs in 
Fig. 1^. The results from the two nondiffusing runs discussed in |^ are also 
included. In all of the diffusing cases, the model overestimates ID by roughly 
the same amount, which we take to be the effect of the overshoots at the 
first two maxima. At the higher temperatures the discrepancy is also higher, 
presumably due to the model's long tail. 

It is interesting to note that this MD system produces results that agree 
very closely with experiment: For example, the MD predicts that at T = 
425.0 K and p = 0.925 g/cm'^, D = 6.40 x 10"^ cm^/s, while experiments by 
Larsson, Roxbergh, and Lodding find that at T = 425.0 K and p = 0.915 
g/cm^, D = 6.020 x 10~^ cm^/s. Hence our agreement with MD results 
genuinely reflects agreement with properties of real liquid Na. 



4 Conclusions 

We have presented a single-atom model of a monatomic liquid that provides 
a unifled account of diffusing and nondiffusing behavior. The nondiffusing 
motion is modeled as a sum of oscillations at the normal mode frequencies 
(Eq. (|T^)), simulating the trajectory of an average atom in a complicated 
single-body potential well that fluctuates due to the motion of its neighbors. 
Self-diffusion is accounted for in terms of instantaneous transits between 
wells, which occur at a temperature-dependent rate u. Since this model 
gives a simple and straightforward account of the motion itself, it can easily 
be used to calculate any single-atom correlation function one wishes; here 
we have focussed on the velocity autocorrelation function. It is interesting 
to note that in this model the velocity correlations persist through a transit, 
instead of being washed out entirely by the transit process; we will return 
to this point below. The relaxation of correlations seen by the decay of 
Z{t) arises here from two distinct processes: Dephasing as a result of the 
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Figure 7: £> as a function of T for both the model and MD. 



large number of frequencies in the single-well motion, and transits between 
wells. The dephasing effect produces relaxation but not diffusion: It causes 
Z{t) to decay [jl| but its integral remains zero. On the other hand, transits 
certainly contribute to relaxation (see 0, where they provide the only re- 
laxation mechanism), but in addition they raise the first minimum of Z{t) 
substantially, increasing its integral and providing a nonzero D. 

The comparison of this model to MD results is generally quite positive, 
particularly for a one-parameter model; the two calculations of Z{t) agree 
strongly over 1000 timesteps. In addition, the match is encouraging over a 
very large range of temperatures, from essentially K to 3Tm- The most 
noticeable discrepancies are the backward drift in the location of the first 
minimum of Z{t), which is present in MD but not the model, and the ten- 
dency of the model to exaggerate certain characteristics of the MD results 
(the maxima at intermediate times and the high-T positive tail). This latter 
effect is responsible for the model's overestimate of the diffusion coefficient, 
though we hasten to add that the model D is still in satisfactory agreement 
with MD results, especially for the liquid at T > (see Fig. |^. 

As in |jl|], it is useful to compare this model and the accompanying results 
to the work others have done using the formalism of Instantaneous Normal 
Modes (INM) and similar methods. Previously, we discussed the advantages 
of our general approach and the superior quality of its results when apphed 
to nondiffusing states; here we will consider matters relating explicitly to 
mechanisms of diffusion. As is noted explicitly by Vallauri and Bermejo 
0, the account of INM by Stratt (see, for example, [^j) does not consider 
diffusion at all; their Z{t) is essentially a sum over cosines, and as such it 
integrates to zero. This is understandable, because as Stratt et al. repeatedly 
emphasize, their approximation is valid only for very short times, so they are 
not attempting to model effects with longer timescales. They compare their 
INM results with MD calculations of Z{t) for states of an LJ system ranging 
from a moderately supercooled liquid to well above the melting temperature, 
and our fits are of roughly the same quality or better in all cases. Authors 
who do attempt to model diffusion usually follow the path suggested initially 



by Zwanzig |jTO|, who thought of the liquid's phase space as divided into 
"cells" in which each atom spends its time before finding a saddle point 
in the potential and jumping from one cell to the next. He imagined as a 
first approximation that the jumps destroy all correlations between the cells; 
since atoms are jumping all the time, he suggested that the net result was 
to multiply the nondiffusing form of Z{t) by a damping factor exp(— t/r) for 
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some timescale r representing the lifetime of a stay in a single cell. Notice 
that Zwanzig provided no dynamical model of the jumping process itself. 
This suggestion has been developed and transformed extensively by Madan, 
Keyes, and Seeley |Tl[], who take a general Zwanzig-like functional form for 
Z[t) and use a combination of heuristic arguments and constraints on its 
moments to specify its dependence on a "hopping rate" u^, the analog of 
for Zwanzig, which they then extract from the unstable lobe of the INM 
spectrum. Although we cannot be entirely sure, as indicated in , we think 
it most likely that their simulations of LJ Ar at 80 K, 120 K, and 150 K 
consist of states comparable to our Figs. ^ and |], and again we would argue 
that our fits are somewhat better. Finally, Cao and Voth [12] also approach 
diffusion by means of a damping factor, and they consider factors of two 
different types, each of which contains parameters that can be determined 
from other calculated or experimental quantities. Their matches with MD 
are actually quite good, but again ours are of at least comparable quality. 

Having claimed that our matches to MD simulations are as good as or 
better than all of the others we have surveyed, let us emphasize a fundamental 
difference between our approaches to diffusion: We provide an account of 
the process of transiting from well to well, so we have a model of the actual 
motion of a mean atom in space that is valid to arbitrary times, and given 
this model we can calculate the effect of transits on correlations. All other 
approaches we know of begin with an account of the motion valid only for 
very short times, calculate Z{t) from this motion, and then try to model the 
effects of diffusion on Z{t) directly, using parameters that are thought to be 
characteristic of jumps between wells. In the diffusing regime, no one else 
we have seen actually calculates {v{t) ■ v{0)) to find Z{t), as we do. In the 
process, we find that some of the assumptions made by others, in particular 
Zwanzig's hypothesis that jumps between wells simply erase correlations, are 
not true. This approach is already yielding insights into the actual motion 
of an atom undergoing a transit. 

Finally, we would like to compare the present mean atom trajectory model 
with an earlier independent atom model by Wallace In developing the 
independent atom model, two arguments were made: (a) the high rate of 
transits in the liquid state shows the need to abandon the normal mode de- 
scription of motion, and instead picture the motion of a single atom among 
a set of fluctuating wells, and (b) the leading approximation to a fluctuat- 
ing well is its smooth time average well. Accordingly, the independent atom 
oscillates with frequency in a smooth isotropic well, and it transits with 
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probability fi to an adjacent identical well at each turning point ||^. What 
we have now learned by considering low temperatures, and especially the 
nondiffusing states of is that the MD system exhibits a velocity decor- 
relation process which results from the presence of many frequencies in the 
single-particle motion, making assumption (b) less reasonable. In the cur- 
rent model these many frequencies are retained, and they are intended to 
represent the strong fluctuations in each single-particle well; in this way the 
current model makes an important improvement over the independent atom 
model, where the well fluctuations were averaged out. Beyond this difference, 
the two models contain similar but not identical treatments of transits, whose 
rate increases with increasing temperature, and which produce self-diffusion. 
The less detailed but simpler independent atom model has proven useful in a 
description of the glass transition ||^ , and the corresponding transit parame- 
ter has been used to relate shear viscosity and self-diffusion in liquid metals 

The next logical step in this work, given the results so far, is to extend the 
mean atom trajectory model and apply it to calculations of more complicated 
correlation functions of the liquid's motion. 
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